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SUMMARY 


Linearized  theory  for  compressible  unsteady  flow  is  used  to  derive 
the  velocity  potential  and  lift  and  moment  coefficients  in  the  form  of 
power  series  in  terms  of  the  frequency  of  oscillation  for  a harmonically 
oscillating  rectangular  wing  moving  at  a constant  supersonic  speed. 

Closed  expressions  for  the  velocity  potential  and  lift  and  moment  coef- 
ficients associated  with  pitching  and  translation  are  given  to  the  sev- 
enth power  of  the  frequency.  These  expressions  extend  the  range  of  use- 
fulness of  NACA  Report  1028  in  which  similar  expressions  were  derived 
to  the  third  power  of  the  frequency  of  oscillation.  For  example,  at  a 
Mach  number  of  10/9  the  expansion  of  the  potential  to  the  third  power 
is  an  accurate  representation  of  the  potential  for  values  of  the  reduced 
frequency  only  up  to  about  0 .08;  whereas  the  expansion  of  the  potential 
to  the  seventh  power  is  an  accurate  representation  for  values  of  the 
reduced  frequency  up  to  about  0.2,  a value  of  this  parameter  large  enough 
to  cover  most  rectangular -wing  flutter  cases  likely  to  occur  at  this  Mach 
number. 

The  section  and  total  lift  and  moment  coefficients  are  discussed 
with  the  aid  of  several  figures . In  addition,  flutter  speeds  obtained 
in  the  Mach  number  range  from  10/9  to  10/6  for  a rectangular  wing  of 
aspect  ratio  4.53  by  using  section  coefficients  derived  on  the  basis  of 
three-dimensional  flow  are  compared  with  flutter  speeds  for  this  wing 
obtained  by  using  coefficients  derived  on  the  basis  of  two-dimensional 
flow.# 


INTRODUCTION 


■# 

A method  for  obtaining  the  air  forces  and  moments  acting  on  harmon- 
ically oscillating  rectangular  wings  in  supersonic  flow  is  given  in  ref- 
erence 1.  This  method  is  based  on  the  expansion  of  the  velocity  potential 
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in  powers  of  the  frequency  of  oscillation  and  is  applied  to  harmonically 
translating  and  pitching  rectangular  wings  to  obtain  expressions  for  the 
associated  forces  and  moments  involving  the  frequency  to  the  third  power . 

In  the  expansion  of  the  velocity  potential,  the  frequency  of  oscil- 
lation and  Mach  number  enter  the  results  in  a combined  form  such  that 
the  range  of  frequency  for  which  the  expansion  to  a given  power  repre- 
sents an  acceptable  approximation  to  the  actual  potential  decreases  as 
the  Mach  number  decreases  toward  unity.  As  a result,  the  force  and 
moment  expressions  of  reference  1 apply  to  a sufficiently  broad  frequency 
range  for  most  flutter  studies  only  if  the  stream  Mach  number  is  greater 
than  about  1.5.  In  order  to  obtain  results  that  cover  a larger  part  of 
the  transonic  range,  that  is,  results  that  apply  over  a wider  frequency 
range  at  Mach  numbers  nearer  unity,  the  present  paper  extends  to  the  sev- 
enth power  of  the  frequency  the  expressions  of  reference  1 for  velocity 
potential,  section  force  and  moment  coefficients,  and  total  force  and 
moment  coefficients.  This  extension  results  in  a coverage  of  frequency 
that  is  generally  sufficient  at  Mach  numbers  as  low  as  about  1.1.  It 
also,  of  course,  increases  the  frequency  range  covered  by  the  approxi- 
mate theory  at  all  supersonic  Mach  numbers . 

Although  the  method  of  reference  1 furnishes  a straightforward 
approach  for  making  this  extension,  the  present  paper  employs  a more 
concise  method  based  on  the  velocity  potential  for  a semi -infinite  wing, 
developed  in  reference  2.  In  this  reference  Stewartson  makes  use  of  the 
Laplace  transformation  to  obtain  the  potential  in  the  form  of  a definite 
integral.  The  integrand  of  this  integral  expanded  in  powers  of  the  fre- 
quency and  integrated  term  by  term  can  be  made  to  yield  results  that  are 
identical  in  form  with  those  obtained  by  the  method  of  reference  1.  For 
the  sake  of  completeness  Stewartson 's  derivation  of  the  velocity  poten- 
tial for  the  semi -infinite  wing  is  reconsidered  herein. 

For  illustration  the  extended  section  coefficients  are  used  to  cal- 
culate flutter  characteristics  for  a rectangular  wing  of  aspect  ratio  4-. 53 
at  several  Mach  numbers  in  the  low  supersonic  speed  range . These  results 
are  compared  with  calculations  made  by  using  two-dimensional  aerodynamic 
coefficients,  obtained  from  reference  3- 


SYMBOLS 


A aspect  ratio 

b one -half  chord 


c 


speed  of  sound  in  undisturbed  medium 
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c2,  section  lift  coefficient 

total  lift  coefficient 

cm  section  pitching-moment  coefficient  about  axis  of  rota- 

tion Xq 

total  pitching-moment  coefficient  about  axis  of  rota- 
tion Xq 

functions  defined  after  equation  (17) 


h 

ho 

• • 

h,a, 

J0(c) 

k 

Li, Mi 
. Li, Mi 
M 

Ma 

Ma 

P 


vertical  displacement  of  axis  of  rotation  xD,  positive 
downward 

amplitude  of  vertical  displacement 

time  derivatives  of  h and  a,  respectively 

Bessel  function  of  zero  order  (first  kind) 

reduced  frequency,  ajb/V 

components  of  section  force  and  moment  coefficients, 
respectively,  defined  in  equations  (29)  to  (32) } 
i = 1,  2,  3,  and.  k 

components  of  total  force  and  moment  coefficients, 
respectively,  defined  in  equations  (37)  and  (38) J 
i = 1,  2,  3,  and  L 

Mach  number,  V/c 

aerodynamic  section  moment  on  wing  about  axis  of  rota- 
tion x0,  positive  leading  edge  up 

total  aerodynamic  moment  on  wing  about  axis  of  rota- 
tion x0,  positive  leading  edge  up 

local  pressure  difference 

aerodynamic  section  normal  force,  positive  downward 


, A1' 


k 
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P 


a 


s 


total  aerodynamic  force  on  wing,  positive  downward 

nondimens ional  radius  of  gyration  of  wing  section  about 
elastic  axis,  ^ Iayrnb^  where  Ia  is  mass  moment  of 

inertia  per  unit  span  about  elastic  axis  and  m is 
mass  of  wing  per  unit  span 

one -half  span  of  wing 


t time 

V velocity 

w(x' ,ym' ,t)  vertical  velocity  at  surface  of  wing  along  chordwise 

section  y'  = ym' 


% 

xo 


x'  = 2bx 
y’  = 2by 
z'  = 2bz 
Zh 
Zm 

z* 

a 

% 

ao 


location  of  center  of  gravity  of  wing  measured  from 
elastic  axis  (see  ref.  3) 

abscissa  of  axis  of  rotation  of  wing  (elastic  axis) 

nondimens ional  rectangular  coordinates  attached  to  wing 
moving  in  negative  x-direction,  referred  to  wing 
chord  2b 


first  bending  mode  shape  of  wing 

vertical  displacement  of  any  chordwise  section  of  wing 
first  torsion  mode  shape  of  wing 

angle  of  attack,  positive  leading  edge  up 

effective  angle  of  attack  due  to  vertical  translation, 

fi/V 


(3  = - 1 


amplitude  of  angle  of  attack  a 
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i = 2y/A 


0h 

phase  angle  between  section  lift  due  to  h 
velocity  h 

and 

9a 

phase  angle  between  section  lift  due  to  a 
position  a 

and 

9hm 

phase  angle  between  section  moment  due  to 
velocity  la 

h and 

9am 

phase  angle  between  section  moment  due  to 
position  a 

a and 

9A 

phase  angle  between  total  lift  due  to  a and 
position  a 

9am 

phase  angle  between  total  moment  due  to  a 
position  a 

and 

K 

density  parameter, 

p 

density  in  undisturbed  medium 

0 

disturbance-velocity  potential 

0) 

frequency  of  oscillation 

% 

first  bending  frequency  of  wing 

% 

first  torsion  frequency  of  wing 

cu  = 2kM2/p2 


ANALYSIS 

Velocity  Potentials  for  Harmonically  Oscillating 
Rectangular  Wings 

» 


As  a first  step  in  the  analysis,  an  integral  expression  is  developed 
for  the  velocity  potential  for  a harmonically  oscillating  semi-infinite 
rectangular  wing  (fig.  1) . As  mentioned  in  the  introduction  this  expres- 
sion Is  given  in  reference  2 and  is  redeveloped  herein  for  the  sake  of 
completeness.  From  the  expression  for  the  semi-infinite  wing,  the  poten- 
tials for  the  various  regions  (see  fig.  2(a))  of  a finite  rectangular 
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wing  are  obtained  either  directly  or  by  appropriate  modifications . The 
analysis  for  the  finite  wing,  like  that  of  reference  1,  is  restricted 
to  the  condition  that  the  Mach  line  from  the  foremost  point  of  one  wing 
tip  does  not  intersect  the  opposite  tip  ahead  of  the  trailing  edge. 

Velocity  potential  for  semi-infinite  rectangular  wing.-  Consider  a 
thin,  flat,  semi -infinite,  rectangular  wing  moving  at  a constant  super- 
sonic speed,  in  a chordwise  direction  normal  to  its  leading  edge  as  shown 
in  figure  1.  The  differential  equation  satisfied  by  the  disturbance- 
velocity  potential  for  the  wing  (when  referred  to  a rectangular  coordi- 
nate system  x',y',z'  moving  uniformly  in  the  negative  x' -direction 
with  the  x'y' -plane  coincident  with  the  mean  position  of  the  wing)  is 

+ v Jl)  0 = (1) 

8*7  ax'2  ay'2  az'2 

where  x'  = 2bx,  y'  = 2by,  z'  = 2bz,  2b  is  the  wing  chord,  and  c 
is  the  speed  of  sound  in  the  undisturbed  medium.  The  boundary  condition 
of  tangential  flow  at  the  surface  of  the  wing,  in  accordance  with  small- 
disturbance  linearized  theory,  can  be  expressed  as 


= v(x',ym’,t) 


dx'  St 


(2) 


where  is  the  vertical  displacement  of  the  ordinates  of  the  surface 
of  any  chordwise  section  of  the  wing  such  as  y'  = ym'  in  figure  2. 

The  wing  is  assumed  to  be  executing  simple  harmonic  motion  with  respect 
to  time  t,  so  that  t enters  only  in  the  exponential  exp(icot),  where 
o>  is  the  frequency  of  oscillation.  Equation  (l)  thus  becomes 


\ 


, d2>|>  = b2  d2!  + 2icoM  b^_  _ o£  . 
by'2  bz'2  bx'2  C Sx'  °2 

where  the  disturbance-velocity  potential  0 is  related  to  \|/  by 


(3) 


0(x',y',z',t)  = \|f(x'  ,y',z')e:^1Jt 


* ^ 

V 

. * 


-V 
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For  the  case  of  the  serai-infinite  wing  having  identical  motion  in  every 
chordwise  section  the  boundary  conditions  that  equation  (3)  must  satisfy 
are : 


' 

X&z' 

~x 

y 

u 

0 

+1 

11 

to 

U V > 0) 

(b) 

itm|f  + V = 0 

ax* 

(y1  < 0) 

(5) 

i|r  = 0 

(x*  < 0) 

(6) 

i)r  — ^ 0 

(z  ' — > ±co) 

(7) 

Equation  (4)  is  derived  from  equation  (2)  and  implies  that  the  normal- 
velocity  distribution  on  the  wing  is  given;  equation  (5)  is  the  condi- 
tion that  the  pressure  be  zero  off  the  tip  of  the  wing;  equation  (6)  is 
the  condition  that  no  disturbances  be  propagated  forward  of  the  wing; 
and  equation  (7)  is  a condition  on  the  behavior  at  infinity  (the  manner 
of  approaching  zero  is  associated  with  the  radiation  condition  of 
Sommerfeld) . Equations  (3)  to  (7)  constitute  the  boundary- value  problem 
for  the  velocity  potential  (f). 

Applying  the  Laplace  transform 


= f e”sx'+(x'  ,y'  ,z')dx' 

J 0 

to  equations  (3)  to  (7)  yields  the  transformed  boundary-value  problem 
in  the  form 


kv  + 1 s* . . ^ (si 

dy'2  dz'2  \ c c7 

= v(s)  (y’  ^ 0)  (9) 

\dz  / z’=±0 


8 


NACA  TN  5076 


\jF  = 0 

(y'  < 0) 

(10) 

0 

t 

1*- 

(z1 — > ±») 

(ID 

The  problem  in.  this  form  is,  as  pointed  out  in  reference  2,  similar  to 
the  problem  treated  in  article  308  of  reference  4.  By  applying  the 
method  of  reference  4 to  the  present  problem,  it  can  be  shown  that  the 
value  for  \|r  at  the  upper  surface  of  the  wing  (z'  = +0)  is  given  by 


For  convenience  equation  (12)  may  be  rewritten  (see,  for  example,  p.  478 
of  ref.  5)  as 


If  use  is  made  of  a table  of  Laplace  transforms  (for  example,  pairs  55 
and  89  of  ref.  6)  and  the  Faltung  or  convolution  theorem,  the  inverse 
of  the  Laplace  transform  in  equation  (13)  may  be  written  as 


where 


^(x',y’,+0) 


X ^ 

if  w(x'  - 6,)G(£,,y,)d6’ 

p Jo 


(14) 


G(x',y') 


U'  < py') 


G(x',y’) 


- p2y,2sec\]^dq 


-i^x’ 

, cp2 


(x'  > py’) 
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and.  = 2b| . Prom  equation  (l4)  the  velocity  potential  in  terms  of 
the  nondimens ional  coordinates  x,  | , and  y becomes 


0(2bx,2by,t)  = - — p— ^ w|2b(x  - |)]G(2b|,2by)d| 


where 


G(2bx,2by)  = J0(fx)( 


(x  < Py) 


G(2bx,2by)  = 


Jofx 


- p^y^sec^  e“iaK 


(x  > py) 


- = 2bai£ 

YP2 

Equation  (15)  is  the  desired  integral  expression  for  the  velocity  poten- 
tial for  a semi-infinite  wing.  Note  that  for  x < py  equation  (15) 
reduces  to  the  potential  for  a two-dimensional  wing  (see  ref.  3) • 

An  expansion  of  the  integrand  of  equation  (15)  in  powers  of  the 
frequency  of  oscillation  is  employed  because  the  indicated  integration 
does  not  appear  to  be  obtainable  in  terms  of  known  functions . The  result 
of  this  expansion  is 


1 q=0  r=0(q  - 2r) 


i — 

r t2\2M; 
)’.[r0  X ' 


Ir(i^y)d| 


Y X 
\ 


i 
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where  r i/si  denotes  the  integral  part  of  q/2  and 


Ir(x,y)  =5  (x  < 3y) 


-1 


Ir(*,y)  = | “, 


sec 


I3y  [1  _ secSjj  dTj  (x  > 3y) 

/ 


It  may  be  shown  by  induction  that  Ir  for  x > 3y  can  be  expressed  in 
terms  of  the  functions  Fn  of  reference  1 as  follows: 


Ir 


> (_i)m  (r  - D* 

x2r  ml(r  - m - 1)  I 


x2(r-m-l)F2(nH.1)(x,y) 


(17) 


where 


Two  relations  involving  Fn  that  are  of  particular  importance  in  the 
next  section  are 


J xmFn(x,y)dx  = — A_(xm+1Fn  - Fm+n+1)  (l8) 

and 


\ 

VE y(x  - iy)  (19) 

where 


Fn(x,y)  = ^ sin“Xy^  + 


n 


Li!  m=l 


ajjjxn-mCpy)®-! 


al  3 


2n 


1 


1 
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2m“1(n  - l)(n  - 2) 
(2n  - l)(2n  - 3)  • 


. (n  - m + 1) 
(2n  - 2m  + l) 


(m>  1) 


Equations  (l6)  to  (19)  provide  the  means  for  obtaining  the  potential  to 
any  desired  power  of  the  frequency  of  oscillation.  They  are  vised  in  the 
next  section  to  develop  expressions  to  the  seventh  power  of  the  frequency 
for  the  velocity  potentials  for  the  various  regions  of  a finite  rectan- 
gular wing  undergoing  torsional  and  vertical  translational  oscillations. 
These  expressions  are  extensions  of  similar  expressions  developed  in  ref- 
erence 1 to  the  third  power  of  the  frequency. 


Velocity  potentials  for  finite  rectangular  wing.-  The  coordinate 
system  and  regions  of  interest  for  the  finite  rectangular  wing  are  shown 
in  figure  2.  On  the  portion  of  the  wing  between  the  Mach  cones  emanating 
from  the  foremost  point  of  each  tip  (region  N in  fig.  2(a))  no  inter- 
action takes  place  between  the  flow  on  the  upper  and  lower  surfaces  of 
the  wing.  On  the  portions  of  the  wing  within  the  tip  Mach  cones 
(regions  T]_,  T2,  and  Tj  in  fig.  2(a))  interaction  does  take  place  between 
the  flow  on  the  upper  and  lower  surfaces.  The  velocity  potential  at  a 
point  in  one  of  these  regions  is  designated  by  0jj,  0^,  0^,  or  0^ 

according  to  the  region  that  contains  the  point. 


For  the  particular  case  of  the  wing  independently  performing  small 
sinusoidal  torsional  oscillations  of  amplitude  Oq  about  some  spanwise 

axis  x0  and  small  sinusoidal  vertical  translations  of  amplitude  hQ, 
the  equation  for  is  (see  fig.  2(b)) 


Zfc  = e3^*  ^(x  - x0)  + hjj  = 2ba(x  - x0)  + h (20) 

Substituting  this  expression  for  into  equation  (2)  gives 


w(2bx,t)  = Va  + 2bd(x  - Xq)  + h (21) 

Equation  (21)  complies  with  the  restriction  inherent  in  the  development 
of  equation  (15) , namely,  that  the  wing  have  identical  motion  in  every 
chordwise  section. 

The  velocity  potential  for  region  Ti  for  the  wing  motion  represented 

by  equation  (20)  is  obtained  directly  from  equation  (l6)  . Putting  the 
expression  for  w(2bx,t)  given  by  equation  (21)  into  equation  (l6), 
•retaining  only  terms  in  the  expansion  up  to  the  seventh  power  of  <5, 
inserting  the  appropriate  values  of  Ir  obtained  from  equation  (17), 
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and.  performing  the  necessary  integrations  with  the  aid  of  equation  (18) 
yields  the  following  form  for  the  velocity  potential  in  region  Tq: 


r~ 

= (h  + Va)Aq(x,y)  + 2bd[x2(x,y)  - x0Aq(x,y)j 


J 


(22) 


where 


M.(x,y)  = | 


Fq  - i EFV,  - ^g(xF2  + 02F5)  + gL  ^C2F2  + (232  - l)F^j 


Jl- 

O) 

kStft 


[x3^32  + 5)f2  - 3xF^  + 23^  - + 7)f2 


96om^ 


30x2F^  + (83^  - hp2  + 3)^5] ^—[3x5  (83^  + 28p2  + 2l) F2  - 


5760Mt 


\ 


10x5  (6p2  + 7^  + 15xF6  + 8p6F, 


iaS7 


11  8064OM6 


7x' 


;(8p^ 


36p2  + 33) F2  - I05x^(2p2  + 3)f^  + 105x2Fg  + 


J 


>>2(x>y)  = f Aq(x,y)dx 
J o 


I 


The  integrations  required  to  obtain  the  function  A2  may  be  readily 
performed  with  the  aid  of  equation  (l8).  The  integrated  values  of  Fn 
(from  eq.  (19))  needed  in  equation  (22)  are  listed  in  appendix  A. 

Use  is  now  made  of  the  potential  equation  (22),  to  obtain  the 

potentials  for  regions  N,  T2,  and  Tj.  The  potential  0^  is  obtained 


I 


i 
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from  0, j by  substituting  y = x/p  in  the  functions  A-j_  and  A2  i*1 

equation  (22),  in  which  case  Fn  becomes 


Jt 

2 n 


and  the  velocity  potential  0^  is  given  by 


(23) 

where 


(h  + Va)f1(x)  + 2bd|f2(x)  - x0fl(^)] 


fn(x)  = X - i5J  4 - — (202  + 3)2#  + ^-(202  + 5)2^  + -^r( 8p4  + 
1 2 ' 5 iad2V  P -7  4 192M^V 


4op2  + 35)^  - -^r-(834  + 56p2  + 63)^ (l6p6  + 1683^  + 

’ 5 96om4  ' ' 6 11520MD ' 

378p2  + 231)^-  + >■  (l6p6  + 2l6p4  + 594  p2  + 429)^ 

7 ' 8o64om° ' ' 8 


f2(x)  = f fq(x)dx 
J 0 

Note  that  because  of  the  two-dimensional  character  of  the  problem  in 
region  N the  spanwise  variable  y is  not  contained  in  0jj.  Equa- 
tion (23)  can  also  be  obtained  by  expanding  to  the  seventh  power  of  55 
the  velocity  potential  for  the  two-dimensional  wing  given  in  reference  3* 
The  potential  0^  is  obtained  from  0ij  by  replacing  y in  the  func- 
tions A^  and  A2  in  equation  (22)  by  A - y,  where  A = 2s/2b  is  the 
aspect  ratio.  The  potential  in  region  T^  (this  region  exists  if 

1 < Ap  < 2)  is  a simple  superposition  of  the  potentials  for  regions  N, 

T]_,  and  T2,  as  indicated  in  reference  1,  and  may  be  written  as 

+ $p2  - 0jj  (24) 


Ik 
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Forces  and  Moments 

Section  forces  and  moments The  expanded  velocity  potentials  for 
the  finite  wing  are  now  used  to  obtain  expressions  for  the  section  forces 
and  moments  at  any  spanwise  station  of  the  wing.  Since  the  distribution 
over  the  entire  wing  is  symmetrical  with  respect  to  the  midspan  section, 
only  expressions  for  the  forces  and  moments  at  any  station  of  the  one- 
half  span  adjacent  to  the  origin  (see  fig.  2)  need  be  considered. 

The  local  pressure  difference  between  the  upper  and  lower  surfaces 
of  the  wing  may  be  written  as 


* - -2p(ir  | + (25) 

The  section  force,  positive  downward,  is  therefore 

P * -2b  f dx  (26) 

^ 0 

and  the  section  moment,  positive  leading  edge  up,  about  the  arbitrary 
axis  of  rotation  x = Xq  is 


Under  the  previously  mentioned  restriction  that  the  Mach  line  from 
one  tip  not  intersect  the  opposite  tip  ahead  of  the  trailing  edge,  equa- 
tions (26)  and  (27)  must  be  evaluated  for  the  two  cases  that  can  arise 
(see  fig.  3)-  These  cases  are:  (l)  The  Mach  lines  from  the  tips  do 

not  intersect  on  the  wing  (A(3  > 2)  and  (2)  the  Mach  lines  intersect  on 
the  wing  but  the  Mach  line  from  one  tip  does  not  intersect  the  opposite 
tip  ahead  of  the  trailing  edge  (l  ^ A(3  ^ 2) . Only  the  final  forms  of 
the  section  forces  and  moments  are  given.  These  forms  are  calculated  by 
deriving  the  pressure  difference  for  the  different  regions  from  the 
appropriate  velocity  potentials,  making  use  of  figure  3 to  determine  the 
limits  of  integration  for  the  regions  involved,  and  using  the  relation 
given  in  equation  (18)  to  perform  the  integrations  involving  the  func- 
tions Fnj  the  results  can  be  written  as 
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P = -4pbV2k2eia>t 


ho 

b 


^1  + ^2)  + + ^4 


= -4pb2V2k^eicu't 


iJa, 


MX  + m2)  + aQ(u3  + iM^ 


(28) 


The  quantities  Lj,  (i  = 1,  2,  3,  and  4)  and  (i  = 1,  2,  3,  and  4) 


are  the  components  of  the  section  force  and  section  moment  coefficients, 
respectively.  The  reduced  frequency  k is  related  to  to  and  55  by 
the  relations 


k = t*±>=  ££-33 
V 2M2 


In  order  to  obtain  self-consistent  expressions  for  the  forces  and  moments, 
the  terms  of  the  velocity  potentials  associated  with  h and  a are 
expanded  to  the  seventh  power  of  53  and  the  terms  associated  with  a 
are  expanded  to  the  sixth  power.  As  a result,  the  coefficients  in  equa- 
tions (28)  are  as  follows: 

Case  l(see  fig.  3(a)):  For  any  section  between  the  tip  and  the 
point  where  the  Mach  line  intersects  the  trailing  edge,  or  where 
0 < y < 1/3, 


In  = 


*33 


4j^32F1  + (202  + i)fJ  - i 3F2  - 63^5  + (804  + 43" 


1 

240 


5(632  + 7)f2  - 30^  - 403%5  + 3(l636  + 804  - 232  + l)f6  - 


20l6o 


7 1634  + W + 33  Fj,  - 35 


.(bp2  + s)f4  + 105F6  - 1123  + 


— 6 

^12838  + 6436  - 1631*-  + 83s  - 5)  Fg 


(29a) 
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L2 


jt0k 


2FX  + 


(232  - l)f2  - 302F3  /g-)2  + i 5F2  - (834  - 402  + 3)^4 


lOp1^ 


,mJ  2880 


21  (23s  + 3)f2  - 7OF4  - 3 (1606  - 
80^  + 602  - 5^f6  + 56p6F7 


"V  1 


(29b) 


Ics  = 


c0M2k2 


2^02  + ljf±  - 0^-L  - (604  + 302  - i)f2  + 02(602  + 5)f^J)2 


^-^02  + 5)  F2  + 1206F5  - (kO06  + 2O04  - 502  + 3)f4  + 


20^(l502  + 13)F5  (g)  - ^ 3(20^  + 2102  + 2l)  F2  - 

io(02  + 7)f4  + 1200^  - 3(ll208  + 5^06  - 1404  + 

702  - 5)f6  + 806(2802  + 25)  F7  (&j6  | - 2XQL! 


(29c) 


l4  = 


jt05k 

702  - 1 


202F1  - (302  + l)F2  + i 3(204  - 02  + l)F2  - 240^F5  + (200^  + 
5(02  + 7)f2  + 10  (806  - + 302  - 3)f4  - 


c5\2  1 


M/  240 
2^0%  + (l6806  + 640^  - 1102  + 3)1 


! . 


if  + sfellK + 27r  + 


33)f2  - 35 (02  + 9)f4  + 2l(l608  - 806  + 6^  - 502  + 5)f6  - 8960^  + 
(5760s  + 23206  - 4401*  + 1502  - 5)ifi  (gf  i - 2xqL2 

\M/  I 


7 

(29a) 


"2  - Jgt  j1*^  - [-(2S2  - i)f2  + (4p2  + l)F^ffij2  + i[l5F2  - 2(8p4  - 

^ + 3)^  ♦ {&£  ♦ S32  - 1)?6](|)lt  - ^^(as2  + i)%  - 105*4  - 

3(i6p6  - 8p^  + 6p2  - 5)f6  + (^P6  + Vkft  - 4p2  + . - 2x^2 

J (30b) 

m3  = ~^2k2  j^P2  + 1)^2  - \ j^P^!  - 3(6p4  + 3p2  - i)f2  + (20P4  + 

21p2  + 3)f^^2  + A ^(p2  + 5)^2  + 32p%  - 2(40p6  + 20p^  - 
5P2  + 3)fk  + (56P6  + 64p^  + lip2  - l)^j(f)4  - ^j^p4  + 

2ip2  + 2l)l2  - 15 (p2  + 7)!^  + l6op^5  - 3(ll2p8  + 56p6  - 

(Equation  continued  on  next  page) 
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1^  + 732  - 5)^5  + (19238  + 23236  + 4431*  - 532  + l)  F8  ^ 
2xq(Mx  + L3  + 2xoLx^ 


(30c) 


= < -8jj32F1  + (232  + 1)  F^j  + t (23^  - 32  + i)f2  - 63^5  + 

0 

32( 532  + ^(ff  - ^(b2  + 7)f2  + - kph  + 332  - - 

900^  + 234(2832  + 19)ftJ(|)4  + J3  (23^  + 2732  + 33)  F2  - 

14  (32  + 9)5^  + 7(163®  - 83®  + 63^  - 532  + 5)f6  - 22k0j  + 

83®  (1532  + ll)  F9  W6  - - 2xo(m2  + L4  + 2xqL2)  (30d) 


where  Fn  (n  = 1,  2,  . . . 9),  given  in  appendix  A,  is  Fn  evaluated 

at  x = 1.  For  any  section  between  the  point  where  the  Mach  line  inter- 
sects the  trailing  edge  and  the  midspan,  or  where  1/3  < y g A/2, 


;;  1_  1 _ ^32  + 3/q)\2  83^  + 2832  + 21/BS\4  _ 

33  ^ 25  \M/  960  W 


6436  + 4323^  + 7923s  + 429 /m\6 
322560  \M/ 


1 + ittL+.7@ 

HM/  192  \M/ 


83^  + 36b2  + 3 

5/s.V 

11520 

\M/ 

(31a) 


(31b) 
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Lx  = 


Y 


3M2k2 


,-.2  j.  ^/— \2  )i  q4  XKq2  ■ — \4 

P 


,2  + 1 _ 0 +.  W\  + 4jf, ,,+  55  p*  4-  35/  <£ 
12  \M/  960  \M, 


836  + l403^  •>•  357P2  + 231/5\6 
8o64o  W 


- 2xqLi 


H = 


35k 


j2  _ ! + P2  ± 5/u>\2  _ 4p4  + 4932  + 63/55\ 
24  \M/  2880  W 


m + 


836  + 1803^  + 56l32  + 429/55^ 
322560  \M/ 


- 2xqL2 


Mx  = 


p3 


4 _ 432  + 5/a5\2  + 83^  + 2832  + 21| 
3 " 15  W 560 


-\4 

a>) 


64p6  + 4323^  + 792p2  + 429 /55\6 
l8l44o  W 


- 2xqLj 


m2  = 


_1_ 

0k 


, 3/£\2  . 5(402  + 7 )lm\k  7(80^  + 3632  + 33)  (m 

W)  "576 fe) 46585 ^ 


-\6 


M*  = 


3M2k2 


j2  + x _ P2  + 3/c5\2  + 434  + 35P2  + 33/feN 

8 Wy  576 


836  + 1403^  + 357 32  + 231(a>'\6 
46080  W 


- 2x0(Mi  + Lj  + 2x0Li) 
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■ it 
t 


(3lc) 


(3ld) 


(32a) 

“ 2xoL2 

(32b) 


(32c) 


r 

■ ' \ 


4 

I v- 

v * 


1 
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+ 32  + 5/g\2  _ 43^  + 493s  + 63/asV* 
15  IM/  " 1680  \M/ 


181440 


(f)  " 2xo(M2  + l4  + ^oLa) 


(32d) 


Case  2 (see  fig.  3(b)):  For  any  section  between  the  tip  y'  = 0 

and  the  point  where  the  Mach  line  from  the  tip  at  y'  = 2s  intersects 

the  trailing  edge  ^or  where  0 < y < A - the  components  of  the  sec- 
tion force  and  moment  coefficients  are  given  by  equations  (29)  and  (30), 
respectively.  For  any  section  between  the  point  where  the  Mach  line 
from  the  tip  at  y'  = 2s  intersects  the  trailing  edge  and  the  midspan 


or  where 


A - ± <y  < 


A/s) 


the  components  of  the  section  force  and 


moment  coefficients,  respectively,  are  obtained  by  first  adding  to  equa- 
tions (29)  and  (30)  the  results  of  substituting  A - y for  y in  these 
equations  and  then  substracting  equations  (31)  and  (32)  from  the  result. 
For  example,  the  term  of  that  does  not  contain  d5  is  given  by 

% = ~ |-432[li(y)  + FX(A  - yj]  + 4(232  + 1)  [l2(y)  + §2(A  - yj|  j-  - ~ 


Total  forces  and  moments.-  Expressions  for  the  total  forces  and 
moments  can  be  derived  by  considering  only  case  1 of  the  previous  sec- 
tion. It  can  be  shown  that  case  2,  a] though  more  cumbersome  to  handle, 
leads  to  the  same  expressions.  Therefore  the  total  force,  positive 
downward,  may  be  written  as 

- . r1/!3  . rA/2 


r vp  rA/2 

P • khJ0  (P,T <P>» 


and  the  total  moment,  positive  leading  edge  up,  about  the  axis  x = xn 
as  0 


fl/3  r A/2 

* khJ 0 4y  + ^ JlU  w„  ^ 
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In  equations  (34)  and  (35)  the  quantities  (P)^  and  (Mq,)^  are  the 

section  force  and  moment,  respectively,  whose  components  are  given  in 
equations  (29)  and  (30),'and  (P)N  and  (^jj  are  the  section  force 

and  moment,  respectively,  whose  components  are  given  in  equations  (31) 
and  (32). 

Upon  performing  the  integrations  indicated  in  equations  (34)  and 
(35)  with  the  aid  of  the  relation 

fl/3l  dv=  g(n  + 2)  1 

J0  n n(n  + 1)  80 
the  results  can  be  written  as 


P = -8pb2V2k2Aeiajt  + iL2)  + c + ilj 

% = -8pb^v2k2Ae  + iM2)  + Oq^Mj  + i%) 

where 
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M.,  = -L.  m2  rOl  + &±- M 


5 


15  m; 


181440 


/ \6 

(m  + Jl-  3 - 

W 3Ap  11 


♦ 


33^  + l6p2  + _ pg  + 10  p1^  4-  24  p2  + i6/d>\6 

240  tMy  4200  \M/ 

2x0^M2  + + 2x^2) 


(38d) 


SOME  PROPERTIES  OF  FORCES  AND  MOMENTS 


Examination  of  some  of  the  properties  of  the  extended  section  and 
total  forces  and  moments  of  the  present  paper  may  be  of  interest.  For 
this  purpose  consider,  first,  the  conventional  coefficients  for  section 
lift  Cj  and  section  pitching  moment  cm  that  may  be  obtained  from 

equations  (28)  as 

=1  - - = k[4k(I.i  + UaK  + *2(l5  + U*)«] 


cm  - - -2  ‘“(Ml  + IM^  + k^M,  + IM^a 


where  = h/v  is  the  angle  of  attack  due  to  vertical  translation. 

From  equations  (39)  the  section  lift-  and  moment-curve  slopes  (complex 
derivatives)  associated  with  vertical  translation  and  pitching  are, 
respectively, 


^ = -i4k(Lx  + iLg) 


E?  - i2k(Mi + "%) 


(40a) 


» 


t 
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and 


dcj 

da 


= 4k2  (Lj  + H4) 


dc, 


m 


da 


= -2k2(Mj  + 1M4)  (40b) 


The  coefficients  c^  and  cm  and  their  associated  slopes  reduce  to 

well-known  results  in  the  steady  case  (k  = 0)  . For  example,  substitu- 
tion of  the  components  given  in  equations  (31)  into  equations  (39)  yields 
the  result  c^  = 4a/{3  (Ackeret's  result)  for  a nonoscillating  two- 

dimensional  wing.  Equations  (39)  and.  (40)  have  the  additional  feature 
that  multiplying  the  components  and  Mj_  (i  = 1 and  2)  associated 

with  vertical  translation  by  k and  the  components  Li  and  Mi 

(i  = 3 and  4)  associated  with  pitching  by  k2  makes  these  quantities 
more  uniform  in  magnitude  and  therefore  more  suitable  for  plotting. 


Now  consider  the  coefficients  for  total  lift  CL  and  total  moment 
Cm  that  may  be  obtained  from  equations  (36)  as 


cL  = - 


2p 


Cm  = = -2 

m 4pbVA 


Jr^r  “ + ^“h  + k2(L3  + 


)^j 


-ik(Mx  + + k2(M5  + iMjj. 


(41) 


From  equations  (4l)  the  total  lift-  and  moment-curve  slopes  (complex 
derivatives)  associated  with  vertical  translation  and  pitching  are, 
respectively, 


* -lH^ + 


dC. 


= i2k(M1  + iMg)  (42a) 


and 


dCj; 

da 


= 4k2(Lj  + iLj 


dC, 


m 


da 


= -2k' 


2(m5  + mk) 


(42b) 


Equations  (40)  may  be  used  to  illustrate  the  extent  to  which  the 
approximate  expressions  of  the  present  paper  for  the  section  force  and 
moment  on  a finite  rectangular  wing  may  be  useful.  A comparison,  based 
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on  these  equations,  is  made  in  figure  4 of  some  exact  and  approximate 
calculations  for  a two-dimensional  wing  vertically  translating  and 
pitching  about  its  midchord  position.  Results  for  a two-dimensional 
wing  are  chosen  because  in  this  case  approximate  and  exact  results  can 
be  compared  and  the  extent  of  convergence  determined.  The  two-dimensional 
case  is  a legitimate  standard  in  the  absence  of  exact  results  for  the 
finite  wing  since  consideration  of  equation  (l6)  will  show  that  this 
equation  is  at  least  as  convergent  when  expanded  to  a certain  power  of 
a>  for  a finite  wing  as  when  expanded  to  the  same  power  of  <5  for  a 
two-dimensional  wing.  In  figure  4 the  components  of  the  section  moment- 
curve  slopes  defined  in  equations  (40)  are  plotted  against  u>  for 
M = 10/9  and  x0  = 0.5>  The  three  dashed-line  curves  appearing  in  each 

of  the  four  parts  of  figure  4 represent  approximate  results  obtained  by 
substituting  the  components  given  in  equations  (52)  into  the  appropriate 
expressions  of  equations  (40)  and  retaining  terms  involving  o>  up  to 
the  third,  fifth,  and  seventh  powers . The  solid-line  curves  represent 
exact  results  obtained  by  using  the  moment  coefficients  tabulated  in 
reference  3*  Curves  of  the  components  of  the  section  lift-curve  slope 
are  not  included  because  they  show  essentially  the  same  convergence 
tendencies . 


As  may  be  noted  in  figure  4,  the  value  of  to  at  which  a particular 
approximation  departs  from  the  exact  theory  is  essentially  the  same 
regardless  of  the  moment  component  considered.  By  extending  the  dP 
results  of  reference  1 ^represented  by  the  curves  in  fig.  k labeled 

"expansion  to  a?,"  except  that  in  the  present  paper  an  term  has 

been  added  to  the  component  2kMq)  to  d£  and  to  <x>J , the  value  of  gd 

at  which  the  approximation  and  the  exact  theory  depart  at  M = 10/9 
has  been  increased  from  0.8  to  1.3  and  2.0,  respectively.  Further  inves- 
tigation will  show  that  these  values  of  to  remain  essentially  the  same 
regardless  of  the  Mach  number  involved.  The  ranges  of  reduced  frequency  k 


|^equal  to 


£ M2  - l\ 
2 M2  j 


in  which  the  various  approximations  adequately  rep- 


resent the  exact  theory  may  therefore  be  expressed  as  follows : 


For  expansion  to  & , 


0 < k < QA.(m2  - l) 


for  expansion  to 


P^l(Mg.,-.„l) 

M2 
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and  for  expansion  to  , 


0 i ki 


From  these  expressions  it  is  readily  seen  that  expanding  to  ffl?  rather 

than  a?  more  than  doubles  the  allowable  k range.  For  example,  at 
M = 10/9  the  limiting  value  of  k is  increased  from  O.O76,  which  is 
too  small,  to  0.19  which  is  sufficiently  large  to  include  most  flutter 
cases  that  would  likely  be  encountered  at  this  Mach  number. 


In  order  to  study  the  effect  of  aspect  ratio  on  the  forces  and 
moments,  consideration  is  given  to  the  spanwise  variation  of  section 
lift  and  moment  coefficients.  The  curves  in  figures  5 and  6 serve  as 
illustrations . In  figure  5 the  magnitudes  of  the  slopes  obtained  by 
substituting  equations  (29)  to  (32)  into  equations  (1*0)  are  plotted 
against  the  spanwise  variable  1 - £ for  the  conditions  M = 1.3, 

A = 4,  Xq  = 0.4-13,  and  u>  = 0.544.  The  quantity  £ is  the  previously 
used  variable  y divided  by  A/2.  The  phase  angles  associated  with 
the  magnitudes  of  the  slopes  shown  in  figure  5 are  plotted  in  figure  6; 
for  example,  the  phase  angle  associated  with  the  magnitude 
dc 


da 


l52  + L42 


is 


0a  = tan” 


1 Ll^Lj.  The  portions  of  the  curves 

in  these  figures  in  the  range  0 1-  £ £ 0 .4  (l/p  < y £ A/2)  also 

apply  to  a two-dimensional  wing  under  the  conditions  listed.  The  effect 
of  aspect  ratio  may  be  noted  in  the  tip  region  0.4  < 1 - £ % 1.0 
(0^  y < 1/p). 


The  phase  angles  0h  and  0^  in  figure  6 are  of  particular  inter- 
est. In  keeping  with  the  fact  that  a component  of  the  force  due  to  h 
in  phase  with  the  velocity  or  a component  of  the  moment  due  to  a 
in  phase  with  the  angular  velocity  d would  be  destabilizing,  values 
of  0^  between  90°  and  270°  and  values  of  0am  between  0°  and  100° 

would  indicate  possible  dynamic  instability  in  pure  bending  and  pure 
torsion,  respectively.  The  coefficient  I*,,  however,  is  always  positive 

so  that  0^  is  between  -90°  and  90°;  this  agrees  with  the  well-known 

fact  that  pure  bending  oscillations  are  always  damped,  at  least  in  poten- 
tial flow.  The  phase  angle  0^*  on  the  other  hand,  is  less  than  180° 

in  the  two-dimensional  region  and  corresponds  to  undamped  or  unstable 
conditions;  whereas  in  most  of  the  region  affected  by  the  wing  tip  the 
phase  angle  0^  is  greater  than  100°,  which  corresponds  to  damped  or 
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stable  conditions.  The  main  feature  of  figure  6,  therefore,  is  that 
aspect  ratio  has  a stabilizing  effect  with  regard  to  the  motion  a. 

This  result  seems  to  indicate  that  a flutter  analysis  which  takes  into 
account  the  spanwise  variation  of  aerodynamic  forces  would  yield  a higher 
flutter-speed  coefficient  than  one  based  on  two-dimensional  aerodynamic 
forces. 

Upon  substituting  equations  (37)  a^d  (38)  into  equations  (42) , the 
overall  effect  of  aspect  ratio  on  the  total  lift-  and  moment-curve  slopes 
can  be  calculated  for  particular  values  of  the  parameters  M,  ffi,  xQ, 

and  A.  Although  the  effect  of  aspect  ratio  A may  change  considerably 
with  small  changes  in  one  or  more  of  the  other  parameters,  some  insight 
into  the  overall  effect  may  be  gained  from  calculations  in  which,  together 
with  A,  one  of  the  parameters  cB,  Xq,  or  M is  varied  while  the  other 

two  remain  fixed.  Figures  7 and  8 show  the  effect  on  dCi^da  and  dCmycLx 
of  varying  A and  ffl  while  keeping  M and  xQ  fixed,  and  figure  9 
shows  the  effect  on  dCmjda  of  varying  A and  Xo  while  keeping  M 
and  5S  fixed.  The  effect  on  the  total  lift-  and  moment-curve  slopes 
of  varying  A and  M can  be  extracted  from  figures  4 and  5 of  refer- 
ence 1. 

In  figure  7 the  magnitudes  and  in  figure  8 the  associated  phase 
angles  of  the  slopes  dC^jda.  and  dC^da  are  plotted  against  as  for 

M = 1.3#  Xq  = 0.5,  and  several  values  of  A.  In  figure  9 the  limiting 
value  as  5>  — > 0 of  the  slope  dC^da  is  plotted  against  Xq  for 

M = 1.3  and  three  values  of  A.  The  curves  in  figures  7,  8,  and  9 
labeled  A = l/p  apply  to  the  least  value  of  A permitted  by  the  fore- 
going analysis  in  that  they  represent  calculations  for  the  combination 
of  aspect  ratio  and  Mach  number  that  causes  the  Mach  line  from  one  wing 
tip  to  intersect  the  opposite  tip  at  the  trailing  edge. 

As  illustrated  in  figure  7(a),  a decrease  in  aspect  ratio  for  a 
given  value  of  co  produces  a decrease  in  the  magnitude  of  the  total 
lift-curve  slope.  This  trend  is  in  general  true  for  the  range  of  as 
considered  and  is  not  affected  by  a change  in  x0  or  M (see  fig.  4 

of  ref.  1).  In  the  case  of  the  magnitude  of  the  total  moment-curve 
slope,  however,  the  trend  with  aspect  ratio  is  greatly  dependent  on  xQ. 

As  may  be  seen  from  figures  7(b)  and  9 (see  also  fig.  5 of  ref.  1),  the 
magnitude  may  increase,  decrease,  or  perhaps  not  vary  at  all  (as  is  the 
case  in  fig.  9 at  Xo  = O.67)  with  a decrease  in  aspect  ratio,  depending 
on  the  value  of  Xq  under  consideration.  This  seemingly  anomalous 
behavior  of  the  total  moment-curve  slope  is  associated  with  the  change 
in  center  of  pressure,  for  prescribed  motions  of  the  wing,  with  change 
in  aspect  ratio.  Also  note  in  figures  7 and  9 that  the  departure  from 
the  two-dimensional  (A  = «•)  results  is  rather  large  for  aspect  ratios 
less  than  about  4. 


i ' 


t _ . 


! 


NACA  TN  3076 


29 


From  figure  9 for  A = 0#  the  following  phase  angles  are  obtained: 
0^  = l80°  for  Xq  < 0.5  and  0^  = 0°  for  Xq  > 0.5.  The  phase 

angle  of  90°  shown  in  figure  8(b)  for  A = »,  5=0,  and  x0  = 0.5 , 
although  associated  with  a slope  of  zero  magnitude  and  not  really  signif- 
icant, is  therefore  a transitional  value.  Note  in  figure  8(b)  that  0^m 

is  less  than  0°  for  A less  than  6;  this  result  also  indicates  that 
aspect  ratio  is  stabilizing. 


APPLICATION  OF  FORCES  AND  MOMENTS  TO  FLUTTER 
Method  of  Flutter  Analysis 


For  the  purpose  of  applying  the  foregoing  results,  a flutter 
analysis  of  the  Rayleigh  type  for  uniform  cantilever  wings  is  now  con- 
sidered. Such  an  analysis  involves  the  selection  of  a set  of  modal  func- 
tions to  approximate  the  flutter  mode,  the  formation  of  the  flutter  deter- 
minant, and  the  solution  of  this  determinant  for  the  flutter  condition. 
(Although  the  use  of  the  force  and  moment  coefficients  derived  herein 
for  the  undistorted  rectangular  wing  results,  in  certain  cases,  in  a big 
improvement  in  accuracy  over  the  use  of  two-dimensional  flow  coefficients, 
perhaps  a more  accurate  but  much  more  cumbersome  procedure  would  be  to 
use  coefficients  derived  for  a distorted  wing.) 

The  flutter  mode  of  the  uniform  cantilever  wing  is  assumed  to  be 
adequately  represented  by  proper  combination  of  the  uncoupled  first 
bending  and  first  torsion  mode  shapes  of  the  wing.  The  flutter  deter- 
minant based  on  these  two  modes  has  the  form 


Ahh 

Aha 

Aah 

Aaa 

m 


where  the  determinant  elements  are  given  by  (see  fig.  10  for  coordinate  £) 


- * ‘ / 1 (k + 


(l*a) 
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v*  « - i 
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*oh  - Vh  dt  - y«r 1(M1  + iM^Vh  d?  (M>c) 

Am  - %2  1 - (irf  A*2«  - j^1  («3  + Ax2«  (**d> 

and.  where  o>  is  the  flutter  frequency.  In  the  event  that  two-dimensional 
air  forces  and  moments  are  used,  the  force  and  moment  coefficients  appear 
as  constants  in  the  integrals  and  can  be  factored  from  under  the  integral 
signs,  and  the  determinant  elements  become 


(it  may  be  noted  that  eqs.  (44)  and  (45)  are  not  restricted  to  any 
particular  Mach  number  range  if  and  (i  = 1,  2,  3>  and  4)  defined 

in  accordance  with  eq.  (28)  are  taken  to  apply  at  either  subsonic  or 
supersonic  speed.)  The  uncoupled  first  bending  mode  shape  and  the 

first  torsion  mode  shape  needed  for  the  evaluation  of  the  aero- 

dynamic integrals  of  equations  (44)  are  plotted  against  1-5  in  fig- 
ure 11.  The  mode-shape  integrals  of  equations  (44)  and  (45)  are  given 


. * A 


NACA  TN  3076 


31 


ZhZ^  d£  = 0.337 


2^2d£  = 0.50 


The  flutter  condition  is  determined  from  the  nontrivial  solution 
of  equation  (1+3)  • This  solution  may  be  obtained  by  various  means  (see 
ch.  XIII  of  ref.  7) . In  the  present  paper  equation  (43)  was  solved  for 
the  unknown  a 3^/%.  For  a particular  wing  and  Mach  number,  for  which 


values  of  M,  tc,  Xq,  Xq,,  ra  , and  are  specified,  the  reduced 

frequency  k was  varied  until  the  ratio  determined  from  equa- 

tion (31)  matched  that  of  the  wing.  In  this  manner  the  values  of  k 
and  ok,/o>  at  flutter  and  consequently  the  flutter-speed  coefficient 

v/boh,  were  determined  for  the  wing  at  the  selected  value  of  M. 


Flutter  Calculations 


The  foregoing  method  of  analysis  was  used  to  obtain  flutter  char- 
acteristics in  the  Mach  number  range  10/9  S M < 20/6  for  a rectangular 

wing  with  the  following  properties:  A = 4.53#  1/k  = 95*3#  xo  = 0.341, 

Xq,  = O.35O,  ra2  = 0.39#  and  = O.583.  Calculations  were  made  on 


the  basis  of  the  components  of  the  section  force  and  moment  coefficients 
given  in  equations  (29)  to  (32).  A numerical  method  for  evaluating  the 

aerodynamic  integrals  resulting  from  the  use  of  these  components  ^for 
example,  J'  + ILg^Z^d^  in  eq.  (44a)  ^ is  outlined  in  appendix  B. 


For  comparison,  calculations  were  also  made  by  using  force  and  moment 
coefficients  for  a two-dimensional  wing  (strip  theory) . In  this  latter 
case  the  determinant  elements  are  of  the  form  given  in  equations  (45) 
wherein  the  components  and  M-j.  are  obtained  from  either  reference  3 


or  equations  (31)  and  (32)  of  the  present  paper.  In  figure  12  the  two 
curves  resulting  from  these  calculations  are  plotted  in  the  form  of 
flutter-speed  coefficient  v/bo^,  against  Mach  number  M.  As  may  be 

noted  in  the  figure,  the  curves  are  well-separated  at  M = 10/9  but 
tend  to  converge  as  the  Mach  number  is  increased.  At  M = 10/9  the 
flutter  speed  obtained  by  using  coefficients  for  a two-dimensional  wing 
is  about  62  percent  of  that  obtained  by  using  coefficients  for  a rectan- 
gular wing,  whereas  at  M = 10/6  it  is  about  95  percent. 
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CONCLUDING  REMARKS 


The  linearized  theory  for  compressible  unsteady  flow  has  been  used 
to  derive  expressions  to  the  seventh  power  of  the  frequency  of  oscilla- 
tion for  the  velocity  potential  and  section  and  total  forces  and  moments 
for  a thin  rectangular  wing  harmonically  oscillating  in  pitch  and  vertical 
translation  and  moving  at  supersonic  speed.  These  expressions  are  exten- 
sions, although  obtained  by  a different  method,  of  similar  expressions 
derived  to  the  third  power  of  the  frequency  of  oscillation  in  NACA 
Report  1028.  As  a result  of  this  extension  the  largest  value  of  reduced 
frequency  for  which  the  theory  is  sufficiently  exact  has  been  increased 

from  k = - — z-ll  to  k = - — =-i,  where  M is  the  Mach  number. 

'M2  M2 

For  example,  at  M = 10 / 9 expansion  to  the  third  power  of  the  frequency 
is  accurate  up  to  about  k = 0.08,  whereas  expansion  to  the  seventh  power 
is  accurate  up  to  k = 0.2.  The  first  of  these  values  for  k is  too 
small  but  the  second  is  probably  large  enough  to  include  most  rectangular- 
wing  flutter  cases  likely  to  occur  at  M = 10/9. 

Flutter  calculations  were  made  in  the  Mach  number  range  from  10/9 
to  10/6  for  a rectangular  wing  of  aspect  ratio  4.53  by  using  the  section 
force  and  moment  coefficients  for  a rectangular  wing  derived  in  the  pres- 
ent paper  and  the  force  and  moment  coefficients  for  a two-dimensional 
wing  given  in  NACA  Report  846.  Comparison  of  the  two  curves  of  flutter 
speed  indicated,  as  may  be  expected,  that  the  use  of  finite-wing  coef- 
ficients is  very  influential  at  Mach  numbers  near  unity  but  becomes  less 
so  as  the  Mach  number  is  increased.  At  M = 10/9  for  the  particular 
wing  analyzed,  the  flutter  speed  obtained  by  using  the  section  coeffi- 
cients for  a two-dimensional  wing  is  about  62  percent  of  that  obtained 
by  using  coefficients  for  a rectangular  wing;  whereas  at  M = 10/6  it 
is  about  95  percent.  Since  flutter  speeds  calculated  on  the  basis  of 
two-dimensional-wing  coefficients  are  generally  very  conservative  in 
the  low  supersonic  speed  range,  it  would  seem  that  flutter  speeds  obtained 
by  using  the  section  coefficients  for  rectangular  wings  of  the  present 
paper  may  compare  more  favorably  with  experiment. 


Langley  Aeronautical  Laboratory, 

National  Advisory  Committee  for  Aeronautics, 
Langley  Field,  Va.,  January  14,  1954. 


NACA  TN  3076 


33 


APPENDIX  A 

SOME  INTEGRATED  VALUES  OF  Fn  AND  Fn 

From  equation  ( 19)  the  values  of  the  functions  Fn  needed  in  equa- 
tion (22)  are  as  follows : 

*1  * |/w(x  - By)  + x Bln'1  /fly/x 

r2  • i/eyCx  - By)  ♦ ~ Bin*1  Vey/x 

y? . &.a  *Jtsa.  g .iB-t 

. ?x?  + + 4 Bln-1^! 

T . l^TTi7)  * 4 Bln-1  ^ 

? 1575  * 

r . ± -rqpyx1*  > 80^.3  96p3y3x2  _♦  isepVx  + j£6fe?  |Ly(;-.  ,yj  + x8 

® 1*158  * 6 

r . 251x6  ♦ 252pyx5.+  2i»|>gyZx4.t.  320SW  + 38l.|)‘>yij»2  ♦ jugg*-*  gStS & l/jly'(x  Tpy)  + *1  .lx."1  \/iyA 
' 21021  ' 7 v 

r . »29i7  + u62«6ay  *■  y»x%V  > 36oxltaV  . aipxVy1*  ■*■  7&e*?tV  * ic2i.»e6y6  + ao^epTyT  ~___  + ^ <ln*ijj^7; 

8 511*03  * B * 

,9 . §aal*j§ggZtt  ±ia£&sg  r...g»»w * aasjg&t  jssa^ad  usages  *jfi8Sgfid  ± g?  . ^ 

The  functions  to  F^  were  given  previously  in  reference  1 but  are 
repeated  here  for  the  sake  of  completeness. 

The  values  of  FQ  needed  in  equations  (29)  and  (30)  are  obtained 

by  substituting  the  value  of  x at  the  trailing  edge  (x  = l)  into  the 
preceding  values  of  Fn.  That  is. 


a 


»n  = Vi.y) 
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APPENDIX  B 

NUMERICAL  METHOD  OF  EVALUATING  AERODYNAMIC 
INTEGRALS  IN  FLUTTER  DETERMINANT 


A numerical  method  of  integration  for  evaluating  the  aerodynamic 
integrals  of  equations  (44)  is  now  outlined.  The  outline  is  based  on 
the  two  cases  discussed  after  equation  (28)  and  depicted  in  figure  3. 

Case  1 (see  fig.  3(a)) • The  spanwise  stations  at  which  values  for 
the  components  of  the  section  force  and  moment  coefficients  given  in 
equations  (29)  to  (32)  are  to  be  found  are  shown  in  figure  10.  Since 
equations  (29)  and  (36)  reduce,  respectively,  to  equations  (31)  and  (32) 

for  £ > 2/Ap,  the  values  of  Fn  needed  to  obtain  the  components  of  the 

section  force  and  moment  coefficients  at  the  various  stations  are  as 
follows : 


V 


n 


Station 

5 

y 

ala~X]/py 

h 

h 

h 

h 

h 

h 

*7 

*6 

?9 

V3y(l  - 3y) 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

1 

*3 

1/43 

Y5a 

*/6 

.95661 

.37005 

.21703 

•15332 

•110J3 

.09640 

.00137 

.07040 

.06204 

2 

2Ap 

1/23 

1/2 

*A 

1.20540 

.55937 

•33958 

.23921 

.18343 

.14046 

.12464 

.10742 

.09440 

3 

34p 

3/43 

VjA 

*/3 

1.40021 

.70402 

.45010 

.32598 

.25329 

.20601 

.17304 

.14562 

.13046 

4,  5,  6 

i4*p 

it/3 

0 

*/2 

1.57000 

.70540 

.52360 

.39270 

.31416 

.26100 

.22440 

.19635 

.17453 

Note  that  these  values  of  Fn  do  not  vary  with  A and  p so  long  as 
Ap  > 2.  From  the  expressions  for  L^  and  Lg  in  equation  (29),  the 
preceding  table  of  values  of  Fn,  the  curve  for  in  figure  11,  and 

the  flutter  parameters 

k = 0.10  A = 4 

P = 0.8307  A0  “ 5Sp  “ °-1505 

x0  = 0.431  A = 0.5  - — = 0.1990 

Ap 


-r-V 


NAQA  TO  3076 


55 


the  following  sample  integration  table  based  on  Simpson's  parabolic 
integrating  rule  may  be  devised: 


Station,  n 

5 

Zh2 

L1  + ^2 

In  = (LX  + iL^Zb2 

0 

0 

1.0000 

0 

0 

1 

.1505 

.6273 

.4549  + 7-l684i 

.2854  + 4.49651 

2 

.3010 

.3469 

•9355  + 9-59321 

.3238  + 3-32811 

5 

.4515 

.1584 

1.3839  + 10.99391 

.2192  + 1.74151 

4 

.6020 

.0520 

1.6646  + 11.62361 

.0865  + .6o42i 

5 

.8010 

.0041 

1.6646  + 11.62361 

.0068  + .04761 

6 

1.0000 

0 

1.6646  + 11.62361 

0 

/ (la  + “2)Zh2«  - Tf(Jb  + Ml  + 2I2  + W5  + III)  + | (H  + Mj  + i6) 

= 0.1456  + 1.66851 

where  ^n  denotes  the  value  of  the  integrand  at  station  n.  In  a 
similar  manner  the  remaining  aerodynamic  integrals  of  equations  (44)  can 
be  evaluated . 

Case  2 (see  fig.  3(b)):  In  this  case,  as  illustrated  by  equation  (33) 

the  components  of  the  force  and  moment  coefficients_for  the  wing  sections 
passing  through  region  Tj  contain  the  functions  Fn(A  - y)  in  addi- 
tion to  the  functions  Fn(y) . As  a result,  there  is  no  arrangement  of 
spanwise  stations  which  will  yield  values  of  Fn  that  do  not  vary  with 

A and  0.  For  this  reason,  it  seems  advisable  simply  to  choose  stations 
at  even  increments  along  the  span  for  the  numerical  integration. 
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(a)  Plan  form  (x'y'- plane). 
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Figure  2.-  Sketch  illustrating  coordinate  system  for  finite  rectangular 
wing  with  two  degrees  of  freedom  a and  h. 
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(b)  Case  2.- Mach  lines  from  tips  intersect  on  wing 
but  not  ahead  of  midchord  ( I $ A>3  $ 2). 


Figure  3»-  Sketch  based,  on  nondimens ional  coordinates  x and  y 
illustrating  different  Mach  line  locations  accounted  for  in 
analysis . 
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(a)  Real  part  of  moment-curve  slope  associated  with  vertical 
translation  of  wing. 


Expansion  to  Gr 

Expansion  to  G)5 

Expansion  to  a)7 


Exact  (ref.  3) 


(b)  Imaginary  part  of  moment-curve  slope  associated  with  vertical 
translation  of  wing. 


Figure  4.-  Comparison  of  two-dimensional  moment-curve  slopes  based  on 
exact  and  approximate  theory  as  a function  of  frequency  parameter  a> 
for  M = 10/9  and  x0  = 0.5. 


Figure  5.-  Concluded 
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Figure  10.-  Sketch  illustrating  coordinates  and  stations  used  in  numerical 

integration  method  of  appendix  B. 


Figure  12.-  Flutter-speed  coefficients  plotted  against  Mach  number  calcu 
lated  from  rectangular-wing  aerodynamic  coefficients  and  from  two- 
dimensional  aerodynamic  coefficients  for  A = 4.53,  1/k  = 95.3, 

x0  = 0.341,  = 0.350,  ra2  = 0.39,  and  ufc/%  = O.583. 
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